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Motivated by experimental results from RHIC, we suggest how a condensate for the 

^ ■ Polyakov loop might produce explosive behavior at the QCD phase transition. This is 

O . due to a rapid rollover of the condensate field below the transition temperature. 
(N ■ 
O ■ 

g '. Many new results have poured forth from RHIC {^/s = 130A GeV). Several appear to be 
^ ! qualitatively different from those at the SPS (i/i = 17A GeV). "Dynamical" fluctuations 
O I in the average transverse momentum, pt, are almost three times larger than at the SPS 
' [|l|. Observed HBT radii are small, ~ 5 — 7 fm 0, and average transverse boost velocities 
are large, ~ .6 c 0. In particular, that the HBT radii Rout/Rside < 1, as a function of 
pion pair pt ■ 200 — > 400 MeV, suggests "explosive" behavior 0. Lastly, suppression of 
r-| ! particles at high pt : 2 6 GeV appears to be a clear diagnostic probe 0J§]. 
j>! I Based upon a condensate model 0, previously we predicted dynamical fluctuations in 
I the average transverse momentum 0] . Using known features of the condensate model [0 , 
;_i ■ here we suggest how it might produce explosive behavior at the QCD phase transition, 
. . and affect the suppression of particles at high pt. 

At infinitely high temperature, by asymptotic freedom QCD is an ideal gas of massless 
quarks and gluons. Lattice data |§] and effective theories p suggest that from infinite 
temperature, on down to perhaps twice the transition temperature, Tc, it really is a 
plasma: a nearly ideal gas of quarks and gluons, augmented with "thermal" masses. 

For temperatures of 2Tc down to Tc, instead of quasiparticles, it is more useful to think of 
a condensate . In a nonabelian gauge theory without quarks, the order parameter for the 
deconfining phase transition is the Polyakov loop, i. This is the trace of the SU{3) color 
Aharanov-Bohm phase factor in the imaginary time direction. By asymptotic freedom, 
the expectation value of i, (i) = io, is one at infinite temperature; from 't Hooft, it is 
zero in the confined phase, T < Tc. Thus in going from 2Tc down to Tc, the condensate 
for the Polyakov loop evaporates. 

In the condensate model, the pressure is determined by a potential for £ |^ , with 

P ~ ^0 Pideal , i^ideal ~ T^ , (l) 



where Pideai is the pressure for an ideal gas of quarks and gluons. In the pure gauge theory, 
where the quarks are quenched, the lattice finds that at 2Tc, p/pideai ~ -8 0, so maybe 
io ~ .95. 
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Lattice data for p/pideai can then be used to fit £o{T). If e is the energy density, 

'-'''T<^4Tif4- (2) 



rp4 Qrp U QJ. 

In the pure gauge theory, the lattice finds that there is a sharp "bump" in (e — 3p)/T^ 
above T^. This bump indicates that the transition occurs in a relatively narrow region of 
temperature. In the condensate model, this refiects a rapid change in the potential for £o- 

Although the lattice data is less reliable with dynamical quarks, it suggests that the 
pressure is "fiavor independent" [Q. In a mean field analysis of the condensate model, 
this implies that the ^-potential is determined mainly by the gluons — for which there is 
good lattice data. Without quarks, the deconfining transition is weakly first order. We 
make the most conservative assumption, that quarks wash out the deconfining transition, 
to leave only crossover. Even so, with dynamical quarks, a bump in (e — 3p)/T^ is still 
evident [^], indicating a rapid change in the £-potential near T^. 

For heavy ion collisions at high energies, the central region exhibits rapid longitudinal 
expansion along the beam direction. For an ideal gas, with boost invariant Bjorken 
expansion the temperature falls as T ~ where r is the proper time. 

If the transition were strongly first order (as occurs for four or more colors), a mixed 
phase exists, and a hydrodynamic analysis indicates that the system lasts for a long time 
at Tc [|l^. Because of this, the HBT ratio Rout/Rside is significantly larger than one. 



2 |rT| . Thinking of the firetube in the central region as a log, this is the "burning log" 



scenario 



11 1 . A strong first order transition implies that there is a mixed phase of bubbles 



at Tc. For slow expansion, bubbles with io = coexist, at equal pressure, with those with 



^0 7^0 



If the transition is crossover, and the potential for i changes slowly, then the system 
evolves smoothly through T^, with a gradual shift of io from ~ 1 to ~ 0. This is a 
"smoldering log", with Rout/Rside > 1- 

If the transition is crossover (or weakly first order), and the potential for i changes 
rapidly, then the system quickly passes through Tc. Suddenly, the system finds itself 
below Tc, at a nonzero value of io, which is no longer the minimum. It then rolls down, 
oscillating about zero, until it settles to = 0. How this roll down occurs is described 
by the following equation: 

^ + ^ - ZsV'i + V{i) + h^<^H = 0. (3) 

Here is the Laplacian in the spatially transverse and rapidity directions, and V{i) is 
the potential for i. The Polyakov loop is coupled to a field $ for Goldstone bosons (vr, 
K, Tj, 1]') through a term ~ h^^'^l [0]. In the condensate model, the potential V{t) is 
completely fixed by the pressure. The coefficient for spatial variations in the Polyakov 
loop, Zs, was taken from that for SU{?>) Wilson lines, Zs ~ 0; actually, Zg ~ g^. 
We also assumed a Lorentz invariant form, Zq = Zg., although probably Zq ^ Zg. 

In we analyzed this equation in the Hartree approximation, where ('^^)) which 

neglects collisions amongst the produced pions. The condensate field £ was assumed to be 
spatially homogeneous, so that only its variation in time mattered. We took 7^ = 0; this 
neglects the effects of longitudinal expansion, which contributes 1/r to 7^, and dissipation. 
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Above Tc, the potential is sharp, and ^ is trapped in a minimum with ^Q ^ 0. For a window 
of temperature which is extremely narrow, AT ^ 2%Tc, the potential changes suddenly 
into that for the symmetric phase: ^ rolls down, and oscillates about zero. The initial 
energy density, which is stored entirely in the I field, is then pumped into the production 
of pions. As the initial energy density is large, the produced pions have large average 
transverse momentum, ~ 5/,^ = 500 MeV. We note that while this agrees with the 
experimentally measured result , in the present model this value cannot be tuned. The 
system is underdamped, and the average field oscillates several times about zero. These 
oscillations produce relatively large fluctuations about the average transverse momentum, 
Pt 0, on the order of ~ 10%. This scenario is like reheating after inflation, or a quench in 
Disoriented Chiral Condensates (DCC's) [0,|13|. It appears unlikely, however, that many 
coherent oscillations of an £-field, sloshing back and forth, generates explosive behavior. 

The STAR results suggest an alternate scenario. Perhaps, due either to rapid longitu- 
dinal expansion, or to an intrinsic term 7^ 7^ 0, one is in a regime where the effects of 
damping cannot be neglected. In this case, the potential evolves quickly, but because of 
the damping, £0 oscillates only a few times. It is conceivable — but not guaranteed — 
that in this case particle production is explosive. Particles are emitted from a shell |]T1|, 

with Rout/Rside < 1- 

A semiclassical calculation of particle production in the condensate model appears 
reasonable. For a large number of colors, Nc, the potential V(£) ~ A*"^ 0, while Zq, Z^, 
7^ ~ 1; thus a semiclassical calculation is nominally valid to ~ ~ 10% for Nc = 3. 



Getting spectra which are exponential in transverse mass is probably easy |T5|; getting 
relative ratios right, such as K/n ~ 0.2, is far from obvious. Explosive behavior should 
also produce a characteristic electromagnetic spectrum, such as in direct photons [|r^ . 

Both scenarios produce relatively large fluctuations about the average transverse mo- 
mentum, which are measurable on an Event-by- Event basis In the underdamped 
scenario, production is from causally connected regions, of size ~ 1 fm |^. Due to Pois- 
son statistics, the fluctuations from each domain are then reduced by the inverse square 
root of the number of domains in one unit of rapidity; If in one unit of rapidity there are 
~ 300 domains, the fluctuations are ~ 10%/\/300 ~ .6%. Further, smaller rapidity bins 
should increase the dynamical fluctuations; for a bin Ay = 0.25, fluctuations increase 
by a factor of two, to ~ 1.2%. This increase in fluctuations on smaller scales is seen in 
the cosmic microwave background. In the damped scenario, if the time scale for damp- 
ing is much less than the spatial correlation length, then fluctuations may be essentially 
constant with the size of the bin. 

CERES at the SPS finds dynamical fluctuations of the mean pt < 3%; STAR finds 
fluctuations ~ 8% Regardless of which scenario applies, the condensate model predicts 
that the fluctuations are concentrated about the mean, and not due to the tails of the 
distribution in pt. Thus the experimental signal should not be significantly affected if one 
cuts on a relatively narrow bin in momentum. STAR included all particles with .1 < pt < 
1.5 GeV; in this model the fluctuations are the same for particles with .3 < pt < .7 GeV. 
Binning in a narrower region in pt also serves to distinguish it from other models. DCC's 
13| or a chiral critical point |jl8[ produce fluctuations at low pt; jets give fluctuations at 



high Pt. 

We can also make qualitative predictions about energy loss p[; in perturbation theory. 
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the energy loss for a fast particle with transverse momentum pt is 

dE/ dx ~ a^moebye/ \ ^ ^ (4) 

where asipt) is the strong coupling constant, moebye is the Debye mass, and is the 
gluon mean free path. In the condensate model, rn^^Yiye ~ ^O' ^"^^ ^o vanishes at Tc. This 
is also seen directly from the lattice data [Q. Assuming that is finite at Tc gives the 
dependence on given above. In QCD, energy loss becomes small as does. 

PHENIX and STAR presented results which appeared to show suppression of particles 
at highpt : 2 — s> 6 GeV [|]. This is seen even by comparing peripheral to central collisions, 
as then no model dependent assumptions need to be made. This suggests that at ^/s = 
130A GeV, one is in a regime in which ^ 0. At the SPS, no such suppression is 
observed, which suggests that £o ~ at ^/s = 17A GeV. (At RHIC, perhaps the largest 
io is still < 1, with £o ~ 1 only at LHC?) 

The outstanding question is then how these various signals change between ^/s = 
17 A GeV and 130A GeV; is it a sudden transition, or gradual 0? If the explosive scenario 
is correct, one might speculate that the transition region is relatively narrow. However, 
the relationship between temperature and ^/s may be far from linear. In saturation models 
1^, for example, T ~ {^/s)'^^'^, so a relatively narrow region in T corresponds to a broad 



region in ^/s. 
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